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Abstract 

We consider chiral perturbation theory in the meson sector at order 
E^. In the terminology of the external field technique, the two-loop 
graphs so generated are of the sunset type. We discuss the evaluation 
of several of these in the case where the masses of the particles running 
in the loops are equal. In particular, we present integral representa- 
tions that are suitable for the evaluation of diagrams in kinematical 
regions where branch points and cuts are present. 



1 Introduction 



In the framework of chiral perturbation theory (CHPT) Green functions 
are expanded in powers of the external momenta and of the hght quark 
masses. The generating functional is constructed by use of an effective la- 
grangian and requires the evaluation of tree graphs at leading order, one-loop 
graphs at next-to-leading order, and two-loop graphs at next-to-next-to- 
leading order. In this article, we describe the evaluation of several two-loop 
graphs in the equal mass case. 

The first complete two-loop calculation in CHPT was performed in Ref. 
[0], in order to investigate the apparent discrepancy of the one-loop predic- 
tion of the cross section 77 — > tt^tt" with the data The topologies of 
the two-loop graphs considered in contain the ones in 

— vector and axialvector two-point functions 

— 71 ^ ewf 

— scalar and vector form factors of the pion (1) 

— 7C7r — > 7171 

— 77 — > TT+TT" 

In other words, knowing how to evaluate the two-loop graphs in 77 — * 7r°7r° 
allows one to calculate those that occur in ([l|) in the equal mass case. There 
are other Green functions where the two-loop graphs have the same topology 
as (1), e.g. TTTT ^ 47r or 7 ^ 47r. The external momenta in those processes 
are, however, in a different kinematical region than in (1), and the integral 
representations worked out below do not apply. 

The evaluation of the two-loop integrals in (1) is not straightforward for 
several reasons: i) CHPT being a low-energy expansion, one has to keep all 
masses at their physical values - the zero mass limit would result in a poor 
approximation of the matrix element, ii) The interaction is of the derivative 
type, which generates polynomials of high degree in the numerator of the 
loop-integrals, iii) In general, one needs the loop-functions in a region where 
branch points and cuts are present. 

Since the work of [0, additional two-loop calculations have been per- 
formed. In the two flavour sector, these are the amplitudes for 77 tt+tt" 
[|], TT ez/7 [P] and TTTT ^ TTTT 0. In the three flavour sector, there exist 
calculations of the vector and axialvector p[ two-point functions, and of 



a combination of vector form factors [lOl. For a review of these calculations 



we refer the reader to 111 
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The calculational methods developed in were apphed in ^, |], ^ |]. 
Further use of them is underway [Q. As these techniques were never made 
pubhc in a coherent manner, we wish to do so here. At the same time, we 
use the opportunity to simphfy the originally used calculational tools. 

We are, of course, aware that this is not the first publication on two-loop 
integrals. Nevertheless, we feel that it would be inappropriate to give an 
overview of what has been previously done in this field, because those cal- 
culations far as we can judge, mostly unrelated to what we aim at 
here. Indeed, in contrast to e.g. the evaluation of two-loop integrals in the 
framework of the Standard Model, where very different mass scales occur, the 
present article deals with applications in SU{2) x SU{2) CHPT, where the 
masses are equal. Furthermore, there is only a limited number of two-loop 
graphs that will ever need to be calculated. We expect that the techniques 
presented below will be useful in this restricted framework, because they rep- 
resent a coherent method to deal with quite different topologies. In addition, 
the same methods can also be applied in chiral SU (3) x SU (3) - where the 
masses are different - see Ref. for the self-energy graph. To give another 
illustration, we expect that the two-loop graphs in K14 decays at order 
can be worked out with these methods in a straightforward (yet admittedly 
tedious) manner. 

The article is organized as follows. In section 2, we elucidate the structure 
of the terms at order in the chiral expansion, in particular the role of the 
two-loop diagrams. The following sections are devoted to the evaluation of 
the self-energy (section 3), the vertex (section 4), the box (section 5) and the 
acnode diagram (section 6). Section 7 contains the summary and concluding 
remarks. The notation is given in appendix A, whereas appendix B contains 
one-loop integrals. The divergences are evaluated and tabulated in appendix 
C. 



2 The diagrams at order E 

The effective lagrangian of QCD in the meson sector consists of a string of 
terms, 

jCejj = C2 + nCi + h^Ce + ■■■ , (2) 

where tree graphs with generate polynomial contributions of order E'^ 
in the energy expansion. These lagrangians contain external sources which 
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Figure 1: Contributions to the generating functional at order E^. The sohd- 
dash hues denote the propagator in the presence of the external fields. Filled 
circles (filled squares) denote vertices from the effective lagrangians C2 i'^i) 
in Eq. (|^). The open square in Fig. e stands for vertices from Cq. Only the 
sunset diagram Fig. a generates genuine two-loop integrals. 



allow one to evaluate the transition amplitudes with the background field 
method. The path integral representation of the generating functional is 




where [dU] denotes the chiral invariant measure. The low-energy represen- 
tation 

Z = Z2 + hZ4 + h^Ze + ■■■ 

is obtained by expanding the lagrangians Cj around the solution of the classi- 
cal equation of motion SJdxC2 = and carrying out the path integral to the 
required order in h. The diagrams which generate the terms of order are 
collected in Zq and displayed in figure |l[ The solid-dash lines stand for the 
propagator in the presence of the external fields. Full circles (full squares) 
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denote vertices from (-^4), whereas the open square stands for a vertex 
from L^. The diagrams at order for a specific process are obtained by 
attaching the external hnes in all possible ways to these graphs. Examples 
are the self-energy, the vertex, the box and the acnode diagram considered 
below (figures @J^J5 and ||, respectively). 

Figure |I|a collects all genuine two-loop diagrams. It is seen that, in the 
language of the external field technique, all two-loop graphs are of the sun- 
set type. The figures |l]b-e display diagrams that amount to products of two 
one-loop integrals, to products of a one-loop integral with a tree graph con- 
tribution from £4, to one-loop graphs with £4, or to tree graphs alone. In the 
following, we reserve the term "two-loop integral" to contributions from the 
sunset graph Fig. |I|a. [There are Green functions where two-loop diagrams 
are completely absent at order - e.g., the vector two-point functions 
The evaluation of these matrix elements then simplifies accordingly.] 

In the following, we outline the evaluation of the two-loop diagrams that 
occur in the process 77 vi^-n^ in the two flavour case, with equal mass for 
the particles running in the loops. 



3 The self— energy 

We evaluate contributions from the self-energy diagram that is displayed in 
figure |[ The case where the masses of the particles running in the loops are 
not identical is discussed e.g. in |TB] with a technique that is very different 
from the one proposed in this work. References to earlier work on the sunset 
graph may be found in [|r^, see also P^|. 
We consider the integrals^ 

(i7;i/^;i/n = (((l;^^;^^^^r) H;^)) , (3) 

with 

D^ = \-l\ , D2 = l-ll , D^ = l-{p-h-l2)\ (4) 
Integration over I2 generates the loop-function 

J{t) = C{w)T{-w) dx[l-tx{l-x)r , t= {p-hY . (5) 
Jo 

^The notation is given in appendix |^. 
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Figure 2: The self-energy diagram. The filled circles denote vertices from 
the effective lagrangian C2 in Eq- (2). The internal lines stand for scalar 
propagators with mass 1. 



The function J{t) is analytic in the complex t— plane, cut along the positive 
real axis for t > 4. We insert the Cauchy representation 



J{t) = r ; -1.5 <w<0 , (6) 

J4 a — t 

and integrate over li. In this manner, we obtain by use of the formulae in 
appendix ^ 

{H-m) = [da]{F2[^2]; (l-x)F2h]pni , 

H^" = 2/4 [do]{2{l-xfF,[z,]p^p^-F,[z,]g^^}^ , 
Z2 = (1 — s(l — x)) X + cr(l — x) ; s = . (7) 
The integration over the variable a in Eq. ([^) converges in the strip 

-1.5 < Re < -1 . 

Using partial integration in x for the last term in H^'^ , the nontrivial integrals 
in Eq. reduce to 

/ [da]{{l-xrF,[z2]}, . (8) 

J4: 
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It remains to extract the finite and infinite parts in Eq. (|^) as w — > 0. 
We subtract and add the first two terms of the Taylor series of F2 [^2] around 
s = 1. The finite part becomes 

/ dal3{{l-xriC2{x,a-s)}, , (9) 

where we have introduced the kernel 

-(15^ {l"fr + (-i)^) . (10) 

It vanishes at s = 1, together with its first derivative, 

lC2{x, cr; 1) = KL^i^x, cr; 1) = , 
and dies off rapidly at large values for cr, 

/C2 = O , cr ^ 00 . 

The infinite part of Eq. (|]) may be expressed in terms of the quantities 

POO 

D{m,n) = / [da]{{l-xrF4y]}, , 

y = + cr(l - a;) , (11) 

that are evaluated and tabulated in appendix |C[ As an illustration of the 
method, we consider the scalar integral H{s). The subtracted function 

R{s) = H{s)-H{l)-{s-l)H'{l) 

stays finite as w ^ 0, 

daP {)C2{x, a -s)}, . (12) 

The poles at w = are contained in H{1) and in H'{1), 
H{1) = D(0,2) 



H'{1) = 2 {/}(!, 3) -D(2, 3)} 

= -C\w)V\-w)^^w + ^w^ + 0{w^)^ . (13) 



The representation (|T2]) is well suited for numerical evaluation at s < 9 only. 
In the region s > 9, H{s) develops an imaginary part. At o? = 4, 



As a result, one has the dispersion relation 

H{-V = / 7 TT^T r , (15) 

that allows one to evaluate H_ also at s > 9. A similar remark applies to the 
Lorentz invariant components of the tensorial integrals and H^'^. 

4 The vertex 

4.1 Tensorial integrals 

Here we consider the vertex diagram Fig. ^ that leads to the integrals 

{{(^■^i'l-^mu^^) ^ (16) 

with 

D, = l-q , D2 = l-iQ- hf , 
Da = 1-/2 , D^=l-{l2 + li-pif , 

Q = P1+P2 , pi = pl = i- (17) 

Integration over I2 gives the loop-function J(t) with i = {pi — li)^, that we 
represent in the dispersive manner (^. We subtract the emerging subdiver- 
gence by writing J{t) = J(0) + Jit). The contribution from J(0) generates a 
nonlocal divergence that is removed by the usual renormalization procedure. 
We do not consider this piece any further and concentrate on the remainder, 

We collect the denominators by J-'[DiD2{a — i)] and rename the Feynman 
parameters for later convenience, 

Xi,X2 X2,X3 , 
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I, Pi 

h 

P2 

Figure 3: The vertex diagram. The filled circles denote vertices from the 
effective lagrangian C2 in Eq. (2). The double external line denotes a current 
of momentum Q, e.g. one scalar or two electromagnetic currents (contact 
term). The internal lines stand for scalar propagators with mass 1. 




such that 



[d<y]\l 



23 



a [\[z,-{h-R)^]\ 
Z3 = a{l - X3) + xly2 , 
7/2 = 1- 53:2(1 - X2) , 

R = {1 - X2)X3Q + {1 - X3)pi , 

s = Q"" = 2piQ = 2p2Q . (19) 

After the shift li li + R the momentum integration may be performed 
with (pl| ), and the tensors V^, V^^ can be expressed in terms of the scalar 
integrals 

Kn[P; = r {P(x2,a;3)F„[^3]}23 , (20) 

where P(x2,X3) is a polynomial in X2,a;3, and where the argument s in 
Vm[P] s] denotes the s-dependence of z^. This procedure automatically gen- 
erates the tensorial structure in the external momenta. The integrals (pO|) 
can be decomposed into the convergent integral V3[P;s] and the divergent 
polynomials Vi[P; 0] , V2[P; 0] by use of the recursion relation 

Kn[P;s] = Vrn[Pi{l,x;)-Q]+mV^+i[s[xl{2x2-l)]Pus] , 
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Piix2,X3) = / Piy,xs)dy , (21) 







obtained from ( pOD by partial integration in X2- Finally, Vi[P] 0] and V2[P] 0] 
may be expressed in terms of the integrals 

E{m,n) = rI^^{(l-x)™F„MK , 

y = + a{l-x) , (22) 

that are evaluated and tabulated in appendix y. In the following subsection 
we illustrate the procedure in case of the scalar integral V{s). 

4.2 The scalar integral 

Performing the above described procedure, the scalar integral V{s) becomes 

V{s) = V3[Ps;s]-iw + 2){EiO,2)-E{l,2)} , 

Ps = xl{y2 + s{w + 2) {1 - 2x2)x2} . (23) 
As w ^ 0, the finite part is 



Vf{s) = lim V,[P,;s]= rf^v{s, a) 
w-*o a 



1 [\l + sx2{l - 3x2))dx2 C ^ . (24) 
m^Y JO Jo z-i 



[IQti^Y Jo Jo 

For 77 — > TTTT , this representation is not well suited, because Vf contains a 
branch point at s = 4, and the physical region for 77 — >■ vrvr is s > 4. This 
branch point manifests itself in a zero in the denominator of the integrand 



in Eq. (p^ along the curve 2:3 = in the square < X2 , xs < 1 . One may 



solve the problem by writing a dispersion relation for Vf. Using 



-=P(-)+m5(^3) (25) 

^3 ^3 



for s — s> s + iO"^, we obtain 



W.,a) = 7j^X^_ dx2X,^ ,.>4 



X2± 



X3+ = ^{a-W^) ,W^ = {a^-4ay2y/' , (26) 
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from where 



1 



dz 



TT J4 Z — S 



■lmv{z, cr) . 



Integration over a gives 



Vf(s) 



1 r°° dz r°° da 



TT J4 z — S Ji 



a 



(3 lmv{z, cr) 



(27) 



The function Vf may be expressed in terms of elementary functions ||T^, 



Vf{s) 



2^2 



TT 



1 

27 



f + —f- 



3sp' 



25 



+ 



TT 



(28) 



with 



/ = P In 



1 + P 



+ in 



1 - 4/s , s > 4 



(29) 



Corresponding expressions hold for any V3[P;s]. In the evaluation of the 
matrix element for e.g. the process 77 ^ tttt , V3 also occurs in the box dia- 
grams considered below. Due to the algebraic complexity of the expressions 
encountered, it may be more useful use representations analogous to Eq. (|27|) . 
They allow for an efficient algebraic treatment. The triple integrals required 
are in any case considerably easier to evaluate than the four-dimensional 
ones used in the box diagrams discussed below. 



5 The box 

We consider integrals of the type 
with 

= l-ll , D2 = l-{k + qif , 
D3 = l-(/i-g2)' , 

= l-ll , D^ = l-{h + h + qi-pif . (31) 

These are generated by the diagram Fig. ||. We consider the case 
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Figure 4: The box diagram. The filled circles denote vertices from the effec- 
tive lagrangian £2 in Eq- (2). The wavy lines stand for the electromagnetic 
current. The internal lines stand for scalar propagators with mass 1. 



?i + ?2 = Pi +P2 , Pi = P2 = 1 , = ?2 = > 

that is relevant for the process 77 — > tttt . Similarly to the vertex diagram 
considered in the previous section, integration over I2 leads us to consider 
the tensors 

t = (Pi-qi-h)' , (32) 

where we again have dropped the nonlocal singularity generated by J(0). 
The parametrization T[D2DiD2,{a — t)] gives 



with 



Z4 — B — Axi , 

A = X2X3{s{l - X2)X3+ {1 -t){l - X3)} = X2X3A , 

B = A + = xl + X2X3{1 - X3){1 - t) + a{l - x^) , 

5 = qiXiX2X3 - q2X3{l - X2) + {qi - pi){l - X3) , 

s = (Pi+P2)',t=(pi-gi)'. (34) 
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The quantity has aheady occurred in the vertex diagram, see (0). With 
the shift 1^1 — 6, the momentum integrations are easily done by use of 
Eq. (prp, and the tensors B^^- '^'^ may be expressed in terms of the scalar 
integrals 

5™[P;s,t] = rI^{P(xi,X2,a:3)F„[z4]}i23 . (35) 

where P(xi, X2, x^) is a polynomial in Xi, X2, X3, and where the arguments s, t 
in -Bm[-P; s, t] denote the s, t dependence of Z4. These integrals are convergent 
at w = for m > 3. We reduce the divergent integrals to the case m = 3 by 
use of the recursion relation 

Bm[P]S,t] = Vm[3X2X3Pi{l,X2,X3)]s]-mBm+l[APi{xi,X2,X3)]S,t] , 

Pi{xi,X2,X3) = / dy P{y,X2,X3) . (36) 





This relation is obtained from (^) by partial integration in Xi. The vertex 
functions Vm have been discussed above, and it remains to determine -83,4. In 
the kinematical region where Z4 7^ 0, these functions may be obtained from 
(|35D via a four-dimensional integration. In the physical region for the process 
77 — i> vrvr , however, vanishes on a two-dimensional surface embedded 
in the hypercube < xi,X2,xs < 1. Analogous singularities occur in the 
physical region for 77r ^ 77r. These zeros in z^ generate branch points at 
s = 4 and at t = 9. We therefore use again a Cauchy representation, and 
consider the region t < 9, where it suffices to use a fixed-t representation 
- the region t > 9 might then e.g. be reached by use of a Mandelstam 
representation. We illustrate the procedure for 

G{s,t) = \imB3[P;s,t]= —(3g{a;s,t), 

1 f P(xi,X2,X3) 1 . . 

By use of (p5|), with z^ — > 2:4, we obtain for the discontinuity of g 

disCs(7(cr; s, t) = (^(cr; s + zO^, t) — (7(0"; s + iO~, t) 
67ri /■^2+ ri xs ^ fB 



(167r 



2^2 



f^^+ f X3 f B \ 
J dx2 J -^P ,a;2,X3j cix3 , (38) 
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Figure 5: The acnode diagram. The filled circles denote vertices from the ef- 
fective lagrangian L2 in Eq. (2). The wavy lines stand for the electromagnetic 
current. The internal lines stand for scalar propagators with mass 1. 



where , x-^j^ are given in (^6]) . Therefore 



G{s,t) = — / —(3disc,g{a;z,t). 39 

zm Ja z — s Ja (J 

In case that the polynomial P contains the variable s, one has to make 
sure to generate the correct asymptotic behaviour through the dispersive 
representation. It may be necessary to pull out factors of s in the numerator 
before doing the dispersive integral - it is in any case useful to check the 
dispersive representation in a region free of cuts by use of Eq. (^). Analogous 
expressions can be obtained for B4, e.g. by first integrating over xi and then 
again using Eq. 



6 The acnode 

We consider the tensorial integral 



^ 1 



A'-' = {^^ir2E^Jj ' (40) 
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that is generated by the acnode diagram Fig. ^, with 

= , D2 = l-{h + qiy , D, = l-ll , 

qi + q2 = P1+P2 , 

ql = ^2=0, p\=pI = 1. (41) 



The first step in the calculation [|T^ is to employ the Feynman parametriza- 
tion as J-'[DiD2] J^[D/iD^] such that 

\ f\ II {k~{l-Xi)qiY{l2-X2q2Y \\ 
2 



A^''' = j dxi I dx 
Jo Jo 



[i-iiy{i-ilY{i-{i,-i2 + kY)l I ' 

k = Xiqi + X2q2 - Pi ■ (42) 

The integrals d'^li can now be performed once the denominators (1 — /^)^ and 
{1 — {h — I2 + kY) are combined. The result is proportional to 

1 [(l-x,Yl2-kY-{l-x,)q'i] 

dX3 X3 



lo ' ' [l-xsil-xsYh-kYV-^^ ■ 
The integrals d'^l2 can be done by combining 

(i-qY{[x,{i~x,)Y'-ii2-kYY"'" 

in the standard manner. The expression A'^'^ then consists of a convergent 
part proportional to tensors built from the external momenta, and a divergent 
piece proportional to g^'^. The convergent part is easy to evaluate numerically 
in the physical region for 77 — >■ vrvr, because it does not contain branch points 
there. In case one wishes to evaluate these terms e.g. at (gi — piY > 9, one 
may again write a dispersion relation for the form factors in question. Here, 
we concentrate on the term proportional to g'^'^, 



A"" = Aog"'' + ■ ■ ■ , 

Ao = -C'{w)Ti-2w)l ['d''x[x,il-xs)]-^x,il-x,y 

2 Jo 



■w ^ 2w 
^5 



(43) 



where 



Z5 = l-X4: + X3(l - X3)xa |i - (1 - x^k"^] ■ (44) 
14 



Expanding the integrand in powers of i/;, we find 

Ao = C{wfr{-wf-{2 + 7w) + —^J^d'xx4lnz5 + 0{w). 



7 Summary and conclusions 



(45) 



1. We have discussed in this article the two-loop diagrams that occur in 
the evaluation of the amplitudes 77 —>■ 7r°7r° in the equal mass case. 
These are the self-energy, the vertex, the box and the acnode graphs, 
displayed in Figs. These two-loop graphs are also the ones that 
occur in the processes listed in Eq. (|l|), and no less. 

2. We first discuss our results for the self-energy, the vertex and the box 
graphs. These diagrams contain, as a subgraph, the one-loop function^ 

that we represent in a dispersive manner, 

[da] 



J{t) 



a-t 



As a result, they may be represented in d dimensions as linear combi- 
nations of the following integrals: 

T{~w~n)j^ [c/a]{Pi2^+"}i Fig. 2 (self-energy) 

r(-^ - n) f {P^zr^"},, Fig. I (vertex) 

Ti-w -n) {^3^r"}^„ Fig. I (box) 

J4 O" J lid 



Here, n denotes an integer, and Pk are polynomials in k Feynman 
parameters. The Zi are polynomials in the external momenta, in the 
Feynman parameters and in a. 



^For the notation, see appendix ^ 
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3. The vertex and box integrals receive further contributions, that contain 
a nonlocal singularity which is generated by a divergent subdiagram. 
These contributions are cancelled by the standard renormalization pro- 
cedure, and we have not considered them further here. 



4. In order to recover the finite and infinite parts in the vertex and in the 
box diagrams at w ^ 0, we have performed partial integrations in the 
Feynman parameters, reducing in this manner the exponent in z'^'^'^. 
The finite parts are obtained by reducing the exponent to n = — 1, 
while the surface terms generated by partial integration produce the 
divergences. In the case of the self-energy diagram, one does not obtain 
a finite result in this manner. We have instead subtracted the first two 
terms of the Taylor series expansion of around s = 1. 

5. In this manner, we are able to express all divergences in terms of the 
integrals 

{D{m,n) ■ E{m,n)} = [da] |l ; {{I - x)^ F^[y]} ^ , 
where 

y = + a{l — x) ] m = 0, 1, 2 . . . ; n = 1, 2, 3 . . . 
We evaluate and tabulate these quantities in appendix y. 

6. In the physical region for 77 tt'^tt^, the vertex and box diagram 
develop branch points and cuts, as a result of which the above rep- 
resentation for the finite part is not appropriate. We instead write 
fixed-t dispersion relations. We provide an integral representation for 
the required absorptive part in each case. 



7. The vertex diagrams can be given in closed form [^, |^. The algebraic 



complexities in the case of 77 tttt suggest, however, that it may 
sometimes be simpler to keep them in the form of the integral repre- 
sentations provided here. Needless to say that this decision is a matter 
of taste. 

8. Finally, we come to the acnode diagram, shown in Fig. ^. Instead of 
presenting the originally used @, |] method, we evaluate it here along 



lines that are similar to the ones suggested in Ref. |T^ for the decay 
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T] TT 77. The kinematics in 77 n n allows for a substantially 



simpler procedure than the one needed in [17 



9. We conclude that, with these methods at hand, one is able to calculate 
many processes at two-loop order in the framework of chiral perturba- 
tion theory. 
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A Notation 

To simplify the notation, we set the pion masses equal to one, 

M^± = M^o = 1 . 
As is customary, we use dimensional regularization and put 

d ^ 

where d denotes the dimension of space-time. Loop integrations are symbol- 
ized by a bracket. 

We combine denominators with 

K...a.]-' = + 

030:3 . . . a;Ar_i(l -X2) + ... + ajv(l - x^v-i)]"^ 
= J^[ai...aN]. (A2) 

Here [dx]N stands for the normalized measure 

N 



[dx]N = mlle[x,(i-x,)K-^dx, , 

J[dx]N = 1 , (A3) 

and 9{x) denotes the step function. We abbreviate multiple Feynman inte- 
grals by 

= fdx{...} , 

Jo 

{•••}23 = 2/ dX2 xsdxsi...} , 
Jo Jo 

{■ ■ ■}i23 = 6 / dxi X2dx2 / xldxsl- . .} . (A4) 
Jo Jo Jo 
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Furthermore, we use the measure 



with 



and 



^ cwr(3/2) . 

^ ^ r(3/2 + w) M ^ ^ 



hm [dcrl = —^—^ da 



B One— loop integrals 

In the text we use the loop-functions 
/ 1 



[z - ilY" 



= Fm[z\ , 

Fm-l[z] , 



Xz-llY"! 2(m-l) 

'_WilL\ ^ !r<r + eyel. 
[z-Zf]-/ 4(m-l)(m-2) '""'^'^J 



They are given by 



F^[z] = z-^'-"^C{w f^"'.\ '"^ , m > 1 . 

i [m) 



In particular, 

(F,; F^; F^; F^) = CH x 



X r(— 1 — w)z ; r(— w) ; 



2z 



We also use 



with 



j(o) = cHr(-«;) 



19 



C The integrals Z)(m, n) and E{m, n) 

Here we consider the integrals 

{D{m,n)- E{m,n)}^ j^[da]\l- f^dxil-x^FM (CI) 



where 



y — -\- cr(l — a;) ; m = 0, 1, 2 . . . , n = 1, 2, 3 . . . 



In particular, we determine the divergent parts in D(ni,n < 3) and in 
E{m,n < 2). By partial integration in x, we obtain the recursion relation 

, r(n-w-2)Q(w + 2-n) 
{3+w+m-n)D{m,n) = ^ 

-n{D{m,n + l)-D{m + 2,n + l)} , 

(C2) 



with 



/■oo 

Q{a) = C{w)r{-w) [da] a" 



= c^Hr(-^)r(-i-^-a)il^. (C3) 

An analogous relation holds for E{m,n) , with Q{w + 2 — n) ^ Q{w + 1 —n). 
One may use these recursion relations to express D{m > l,n < 3) and 
£■(777., n <2) through the divergent quantities Q and the convergent integrals 
D{m > 1,4) and E{m,S). The case D(0, n) must be treated separately, see 
below. 



C.l Explicit expressions for D{m,n) 
Let 

D(m, n) =C^(u')r^(— w) |p(m, n, 0)+u'p(m, n, l)+w^p(m, n, 2) + 0(w^)| . 

(C4) 

For m> 1, we proceed as described above and find 

p(m, 1, 0) = (m^ + 4m + 5)m(m + 4)A^4, 
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p{m, 1,1) = -(2m^ + 29m^ + 172m^ + 540m^ + 964m^ 
+951m^ + 430m + 36)m(m + 4)A/'2, 

p(m, 1,2) = {m(m + l)D(m + 6) — 3m(m + 3)D(m + 4) 

+ 3(m + l)(m + A)D{m + 2) - (m + 3)(m + 4)D(m)} A^4 
+ (3m^2 + 74m" + 812m^° + 5230m^ + 21938m^ 
+ 62724m^ + 123986m^ + 167682m^ + 149409m^ + 81146m3 
+ 23372m2 + 3456m + 864) m(m + A)N^, 

p(m, 2,0) = -m(m + 2)A^2, 

p(m, 2, 1) = (2m^ + 7m^ + 7m + l)m(m + 2)A^2 ? 

p(m,2,2) = {mD(m + 4)-2(m + l)D(m + 2) + (m + 2)L)(m)}A/'2 

- (4m^ + 26m^ + 61m^ + 63m^ + 27m2 + 4m + 2) m(m + 2)NI 

p(m,3,0) = 0, 

p(m,3, 1) = — l/(4m), 

p(m,3,2) = {D(m + 2) -D(m)}/(2m) - (2m- l)/(4m2). 



(C5) 



where 



N^^ = m{m+l){m + 2), N^^ = N2^im + 3){m + A), 

rco rl dx{l - x)™ 



POO r 

D{m) = / dap 

J4 JO 



{x2 + a(l-x)}^ 
For example, 

D{1) = (tt^ -4)/16, D(2) = -(tt^ - 12)/16, D(3) = -(1371^ - 132)/48 . 

(C7) 

In table 1, we display some of the coefficients p{m, n, k) for convenience. We 
now turn to D{0, n) which is divergent for any n. The recursion relation (|C^ ) 
allows one to evaluate D[0, 1) and D{0, 3) from 

D{0,2) =C(w)T{-w) [da] J dx [x^ + a{l - x)) , (C8) 
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and from D{m > 2,n). In order to evaluate D{0,2), we add and subtract 
from the integrand in ( |C8|) the quantity 

A = (x + a(l - x)r - — - x)r . 

a 

The integral 

C{w)T{-w) J^[da] 1^ dx ((x^ + a{l - x))"" - A) (C9) 
is finite at w = 0, whereas the divergence is contained in 

C{w)T{-w) [da] dxA. (CIO) 
Ji Jo 

In this manner, we obtain the values p(0, n < 3,k) displayed in table 1. 

C.2 Explicit expressions for £'(m, n) 
Let 

E{m, n) = C^(w)r^(— w) |g(m, n, 0)+wg(m, n, l)+w^g(m, n, 2) + 0(w'^)| . 

(Cll) 

Proceeding in the manner described above, we find 
g(m, 1,0) = (m^ + 4m + 2)A^3, 

g(m, 1,1) = -(3m^ + 31m^ + 124m^ + 235m2 + 205m + 64)A/'|, 
g(m, 1,2) = {{m + l)E{m + A)-2{m + 2)E{m + 2) + {m + ?,)E{m)]N^ 
+ (jm^ + 115m^ + 802m^ + 3097m^ + 7230m^ 

+ 10425m^ + 9041m2 + 4295m + 848) NI, 
g(m,2,0) = iVi/2, 
g(m,2,l) = (2m + l)Ar2/2, 

g(m,2,2) = {E(m + 2) - E(m)} iVi + (4m^ + 6m + 3)A^i /2, 

(C12) 

where 

iVfi = (m + 1), A/'g-^ = A^f^(m + 2)(m + 3), 
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E{m) ^ —p 

J A a Jo 



For example 



da /■! dx{l — xy 



x"^ + (t(1 — x) 



1 TT^ — 9 47r2 — 39 



For convenience, we display some of the coefficients q{m, n, k) in table 2. 
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Table 1: The coefficients p{m, n, k). 



m n p{m,n,0) p{m,n,l) p{m,n,2) 






1 


1 

J-O 

12 


144 


1728 


1 


1 


5 
12 


781 
720 


78121 
43200 


2 


1 


17 
60 


2389 
3600 


233857 
216000 





2 


3 
2 


17 
4 


59 
8 


1 


2 


1 
2 


17 
12 


59 
24 


2 


2 


1 
3 


59 
72 


1333 
864 





3 


1 
4 


1 
2 


6 ^ 


1 


3 





1 
4 


7r2 1 5 
6 4 


2 


3 





1 
8 


7r2 1 23 
6 16 


3 


3 





1 

12 


7r2 , 3 
6 2 


4 


3 





1 

16 


7r2 1 883 
6 ~'~ 576 
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Table 2: The coefficients q{m,n,k). 



m n q{m,n,0) q{m,n,l) q{m,n,2) 
1 



1 1 

2 1 

2 

1 2 

2 2 

3 2 

4 2 



1 _ 16 223 

3 9 54 

7_ _331 9011 

24 288 3456 

7_ 1499 206087 

30 1800 108000 

1 1 27if _ n 

2 2 3 2 

1 3 271^^ _ 93 

4 8 3 16 

1 _5_ 2tt^ _ 325 

6 18 3 54 

1 7_ 22ll — 7073 

8 32 3 1152 

J_ A 27!f _ 27983 

10 50 3 4500 
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